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Abstract 

A highly efficient numerical algorithm for solving 3-D 
hypersonic flows with applied electromagnetic fields is 
being developed. Maxwell’s equations of electromag- 
netodynamics have been added to the PNS equations 
and the resulting set of equations are solved using a 
loosely-coupled approach. The fluids and magnetic 
field equations are decoupled in such a way that the 
physics of the problem is retained. To account for up- 
stream (elliptic) effects, the flowfields are computed 
using multiple sweeps with the iterated PNS (IPNS) 
algorithm. The fluid flow and magnetic field solvers 
can either be updated independently or coupled iter- 
atively which makes for an efficient approach in deal- 
ing with the effects of the magnetic field on the fluid 
flow. The new algorithm has been used to determine 
the influence of magnetic fields on two-dimensional 
hypersonic flows over flat plates and wedges. The 
present results are in good agreement with previous 
Navier-Stokes calculations. 

Introduction 

The recent renewed interest in magnetohydrodynam- 
ics (MHD) can be attributed to at least two potential 
applications. The first is the possibility of reducing 
supersonic drag and heat transfer by imposing a mag- 
netic field on the ions produced at the bow shock 
wave [1]. The second is the possibility of increas- 
ing the efficiency of hypersonic air-breathing engines 
by using the energy bypass concept [2]. The MHD 
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energy bypass concept is based on redistributing en- 
ergy between the various stages of a ramjet/scramjet 
engine in order to reduce the Mach number at the 
entrance to the combustion chamber. 

Flowfields involving MHD effects have typically 
been computed [3-8] by solving the complete Navier- 
Stokes (N-S) equations for fluid flow in conjunc- 
tion with Maxwell’s equations of electromagnetody- 
namics. When chemical, vibrational and electronic 
nonequilibrium effects are also included, the compu- 
tational effort required to solve the resulting coupled 
system of partial differential equations is extremely 
formidable. One possible remedy to this problem is 
to use the parabolized Navier-Stokes (PNS) equations 
in place of the N-S equations. The PNS equations can 
be used to predict three-dimensional, hypersonic vis- 
cous flowfields in a very efficient manner [9]. This 
efficiency is achieved because the equations can be 
solved using a space-marching technique as opposed 
to the time-marching technique that is normally em- 
ployed for the complete N-S equations. 

One of the most widely-used PNS codes is NASA’s 
upwind PNS (UPS) code which was originally devel- 
oped by Lawrence et. al [10]. The UPS code solves 
the PNS equations using a fully conservative, finite- 
volume approach in a general nonorthogonal coor- 
dinate system. The UPS code has been extended 
to permit the computation of flowfields with strong 
upstream influences. In regions where strong up- 
stream influences are present, the governing equa- 
tions are solved using multiple sweeps. As a result 
of this approach, a complete flowfield can be com- 
puted more efficiently (in terms of computer time 
and storage) than with a standard N-S solver which 
marches the entire solution in time. Three itera- 
tive PNS algorithms (IPNS, TIPNS, and FBIPNS) 
have been developed. The iterated PNS (IPNS) al- 
gorithm [11] can be applied to flows with moder- 
ate upstream influences and small streamwise sepa- 
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rated regions. The time iterated PNS (TIPNS) algo- 
rithm [12] can be used to compute flows with strong 
upstream influences including large streamwise sep- 
arated regions. The forward-backward sweeping it- 
erative PNS (FBIPNS) algorithm [13] was recently 
developed to reduce the number of sweeps required 
for convergence. 

In the present study, the UPS code has 
been extended to permit the computation of 
two-dimensional/axisymmetric MHD flowfields. 
Maxwell’s equations of electromagnetodynamics 
have beeen added to the PNS equations and the 
resulting set of equations are solved using a loosely- 
coupled approach, similar to that used previously 
to compute flows in chemical, vibrational, and elec- 
tronic nonequilibrium [14,15]. To enhance coupling 
and also to account for the elliptic character of 
the MHD equations, the flowfields are computed 
iteratively using multiple streamwise sweeps with the 
IPNS algorithm. The resulting code has been tested 
by computing two cases involving fluid flows with 
applied electromagnetic fields. These cases include 
viscous flows over flat plates and inviscid flows over 
wedges. Comparisons have been made with the 
previous N-S computations of Hoffmann et al. [7] 


Governing Equations 

The governing equations for a viscous MHD flow 
are the Navier-Stokes equations for the flow of 
the fluid, along with Maxwell’s equation of elec- 
tromagnetodynamics. The equations are given by [7]: 

Continuity equation 

^ + V • (pV) = 0 (1) 


Magnetic field conservation equation 


V • B = 0 (5) 

The flow is assumed to be in chemical equilibrium, 
and the curve fits of Srinivasan et. al [16, 17] are used 
for the thermodynamic and transport properties of 
equilibrium air. Powell’s source term [7] is added to 
the governing equations in order to eliminate the sin- 
gularity due to the existence of zero eigenvalues of the 
flow system. The source term, which is proportional 
to the divergence of the magnetic field is given by 


H = H M (V-B) (6) 

where H^v/ is defined later. Since the divergence of 
the magnetic field is zero, the addition of this source 
term does not change the governing equations. 

The governing equations are nondimensionalized 
using the following reference variables. 




L 


* * 

U - — } t 

U co 


Uoot 

L 


P* = 


P 

Poo 


e, = 


et 

U 2 

u oo 



PooUc* Poo 


n* D* _ & x, By 

z} y t/oov/AW^ 


(7) 


Pm = 


Pm 
Pm co 


= 1 


= 


Vm c 


Vm c 


where the superscript * refers to the nondimensional 
quantities. In subsequent sections, the asterisks are 
dropped. The magnetic Reynolds number is defined 
as 
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Momentum equation 
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Energy equation 


The magnetic field equations are then uncoupled from 
the governing equations of the flow, resulting in one 
set of equations for the fluid flow and another set of 
^ ^ governing equations for the magnetic field . 

( 2 ) 

Fluid Flow Equations 
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The, vector form of the nondimensionalized governing 
equations for fluid flow (with magnetic field effects) 
in a 2-D Cartesian coodinate system are given by 


(3) 


du , aE,- . <9 f,- „ _ dE„ , dF„ 

dt dx dy dx dy 


( 8 ) 


Maxwell’s equation 


The vector of the dependent variables is 


- q £ + V x (V x B) = t^ m V 2 B (4) 


T 

u = [p,pu,pv,pe t \ 


(9) 


2 

American Institute of Aeronautics and Astronautics 


where 
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and the in viscid flux vectors are given by 
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Powell’s source term for the fluid flow is given by 
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where the nondimensional shear stresses and heat 
fluxes are defined in the usual manner [7,9]. 


and the inviscid and viscous flux vectors are given by 
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Powell's source term for the magnetic field is 
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Generalized Coordinates 

The governing equations for the fluid flow and the 
magnetic field are transformed into computational 
space and written in a generalized coordinate system 
($>*?) as 

iu { + E f + F„ + j = 0 (25) 

1 

-B, + Ef + F™ + — = 0 (26) 

where 

E = (E» — E v ) + (F» — F„) 

F = (y-j (E< — E„) + (y-) (F - F„) (27) 

E m = ^ (Er - e?) + (Fr - k) 

F m = (^) (E, m - E") + (&) (Ff - F") 

PNS Equations 


The fluid flow equations with magnetic field effects 
Magnetic Field Equations [Eq. ( 25 )] are parabolized by dropping the time 

The vector form of the nondimensionalized governing derivative term and the streamwise (£ direction vis- 
equations for the magnetic field in a 2-D Cartesian ecus fkw terms in the flux vectors. Equation (25) 
coordinate system are given by 

dB <m? dF? d ei dF£ 

dt dx dy dx dy 
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The solution vector for the magnetic field is 
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The prime in the preceding equation indicates that 
the stream wise viscous flow terms have been dropped. 


Numerical Solution 

The MHD flow is solved using a loosely-coupled ap- 
proach, wherein the fluid flow and magnetic field 
problems are separated from each other. The re- 
sulting system of equations are coupled through the 
momentum and energy equations. The effects of the 
magnetic field are manifested in the momentum equa- 
tions through the electromagnetic forces, and in the 
energy equation. The effects of the velocity field ap- 
pear in the convective terms of the magnetic field 
equations. 

The IPNS flow solver is first used to obtain a veloc- 
ity profile at each streamwise station. This solution 
is then used in the magnetic field equations to deter- 
mine the magnetic field profile. The above process is 
repeated at each streamwise station until a complete 
(initial) sweep of the flowfield is made. Subsequent 
sweeps of the flowfield utilize downstream informa- 
tion from the previous sweep (at each station) to ac- 
count for upstream (elliptic) effects. This iterative 
procedure is continued until a converged solution is 
obtained. 


Fluid Flow Solver 

The iterative PNS (IPNS) method [11] has been mod- 
ified to solve the PNS equations with magnetic field 
terms included, i.e. Eq. (28). The E vector is split 
using a modified Vigneron parameter w m . This pa- 
rameter accounts for the magnetic field effects on the 
ellipticity of the flow, and is given by 


uj m ~ min 
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7 'Ml 
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(30) 


where M x = uj a and a is the nondimensional speed 
of sound. This parameter is derived assuming B x = 0 
and v << u. As a result, E can be written 
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The streamwise derivative of E is then differenced 
using a forward difference for the “elliptic” portion 

(E*): 

(^),. +1 = h. l(E: « ‘ En + <Ef « ■ E? +‘» 

(33) 

where the subscript (i + 1) denotes the spatial index 
(in the £ direction) where the solution is currently 
being computed. The vectors E* +1 and Ef +1 are then 
linearized in the following manner: 
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where 


E = E* + E p (31) After substituting the above linearizations into 

Eq. (33), the expression for the streamwise gradient 
of E becomes 
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The final discretized form of the fluid flow equa- 
tions with magnetic field effects is obtained by sub- 
stituting Eq. (36) into Eq. (28) along with the lin- 
earized expression for the flux in the cross flow plane. 
The final expression becomes: 





- fi+ 1 

(U; + i - U,) fc + 1 


= RHS 


(37) 


where 

RHS 
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Test Case 1: Hypersonic flow over a flat 
plate with applied magnetic field 

This test case corresponds to the flat plate test case 
of Hoffmann et al. [7], wherein the air flow is hyper- 
sonic, viscous, and in chemical equilibrium. The flat 
plate is maintained at a constant temperature and a 
fixed magnetic field is applied at the inflow boundary. 
Along the flat plate, the magnetic field is determined 
by specifying that the normal gradients of the mag- 
netic field components are zero. Likewise, zero gra- 
dients of the magnetic field components are assumed 
at the top and outflow boundaries. A schematic of 
this test case is given in Figure 1. The dimensional 
flow parameters for this test case are 
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drj 
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and the superscript k+1 denotes the current iteration 
(i.e. sweep) level. 
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Magnetic Field Solver 

The magnetic field equations (Eq. (26)] are solved 
at each streamwise station (i 4- 1) after the fluid 
flow equations have been solved using the algorithm 
described in the previous section. Thus, the ve- 
locity components, u and v ) are known at station 
(i + 1). The unsteady magnetic field equations are 
then solved using the Euler implicit scheme [9]. Typ- 
ically, 5-10 time steps are required to converge the 
solution at each station. This solution yields the B x 
and By profiles in the crossflow plane at station (i+1). 
These values of B x and B y can then be put back into 
the fluid flow solver to recompute values of u and v 
at the current station in a predictor-corrector fash- 
ion. These corrected values of u and v can then be 
used in the magnetic field solver to recompute values 
of B x and B y at the ( i -1- 1) station. 

Numerical Results 

In order to investigate the utility and accuracy of the 
present approach of solving an MHD flowfield in a 
loosely-coupled fashion, a few basic test cases were 
computed. The hypersonic flow in these cases was 
altered by the presence of the magnetic field which is 
applied to the flow. 


For this case, an extremely large constant value of 
electrical conductivity is specified which makes the 
fluid nearly infinitely conducting. It should be noted 
that Hoffmann et al. [7] erroneously give a value for 
the electrical conductivity of 100 mho/m for this test 
case. 

A highly stretched grid consisting of 100 points in 
the normal direction was used to compute this case. 
The first point off the wall was located at 5 x 10“ 5 * 
m. Figure 2 shows a comparison of the x-component 
velocity profile at x/L = 0.99 for different values of 
Bq . The present results compare very well with the 
Navier-Stokes results of Hoffmann et al. As can be 
seen in the figure, the boundary layer grows as the 
magnetic field strength is increased. A comparison 
of the temperature profiles at x/L = 0.99 is shown in 
Fig. 3. The present results agree closely with those of 
Hoffmann et al., although slight differences can be de- 
tected even for the case of no applied magnetic field. 
The present results were computed using equilibrium 
air curve fits [16, 17] for both the thermodynamic and 
transport properties and this may explain the small 
differences. Profiles of the magnetic field components 
B x and B y are compared in Figs. 4 and 5 at x/L = 
0.99. The present results are in reasonable agreement 
with the N-S results of Hoffmann et al. 
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Test Case 2: Hypersonic flow over a 

wedge with applied magnetic field 

The second test case considered involves the steady 
hypersonic 2-D flow of an inviscid, infinitely conduct- 
ing fluid past a non-conducting wedge. This test case 
corresponds to the one solved previously by Hoffmann 
et al. [7] A rectilinear, uniform magnetic field which is 
constant along the surface of the wedge is applied in a 
direction orthogonal to the freestream. A schematic 
of this test case is given in Fig. 6. The flow parame- 
ters for this test case are 
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= 5° 

M c o 

= 10.0 

Poo 

= 299.8 N/m 2 

Too 

= 260.91 K 

[B x ) wall 

= 0.0 

(By) wall 

= B 0 = 0.05 T 


The pressure and temperature profiles are pre- 
sented in Figures 7 and 8, for zero angle of attack. 
Comparisons are made with the numerical solution 
of Hoffmann et al. [7]. It can be seen that there is 
good agreement between the current pressure results 
and the results of Hoffmann et al. Computations in- 
volving larger values of Bo are currently underway. 

Concluding Remarks 

Although only limited results have been obtained 
thus far, it can be seen that the present approach 
of decoupling the magnetic problem from the fluid 
flow problem, and solving them in an iterative man- 
ner is quite promising. Computations of other test 
cases are currently underway in order to validate the 
accuracy of the current method in a variety of flow 
situations including internal flows. 
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Figure 1: Test Case 1 



Figure 2: Velocity profiles 
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Figure 3: Temperature profiles 
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Figure 4: Magnetic field (x-component) profiles 
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Figure 5: Magnetic field (y-component) profiles 
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Abstract 

A new parabolized Navier-Stokes (PNS) algorithm 
has been developed to efficiently compute magne- 
tohydrodynamic (MHD) flows in the low magnetic 
Reynolds number regime. In this regime, the electri- 
cal conductivity is low and the induced magnetic field 
is negligible compared to the applied magnetic field. 
The MHD effects are modeled by introducing source 
terms into the PNS equation which can then be solved 
in a very efficient manner. To account for upstream 
(elliptic) effects, the flowfields are computed using 
multiple streamwise sweeps with an iterated PNS al- 
gorithm. Turbulence has been included by modify- 
ing the Baldwin-Lomax turbulence model to account 
for MHD effects. The new algorithm has been used 
to compute both laminar and turbulent, supersonic, 
MHD flows over flat plates and supersonic viscous 
flows in a rectangular MHD accelerator. The present 
results are in excellent agreement with previous com- 
plete Navier-Stokes calculations. 

Introduction 

Flowfields involving MHD effects have typically been 
computed [1-10] by solving the complete Navier- 
Stokes (N-S) equations for fluid flow in conjunction 
with Maxwell’s equations of electromagnetodynam- 
ics. When chemistry and turbulence effects are also 
included, the computational effort required to solve 
the resulting coupled system of partial differential 
equations is extremely formidable. One possible rem- 
edy to this problem is to use the parabolized Navier- 
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Stokes (PNS) equations in place of the N-S equa- 
tions. The PNS equations can be used to compute 
three-dimensional, supersonic viscous flowfields in a 
very efficient manner [11]. This efficiency is achieved 
because the equations can be solved using a space- 
marching technique as opposed to the time-marching 
technique that is normally employed for the complete 
N-S equations. 

Recently, the present authors have developed a 
PNS code to solve supersonic MHD flowfields in 
the high magnetic Reynolds number regime [12]. 
This code is based on NASA’s upwind PNS (UPS) 
code which was originally developed by Lawrence et 
al. [13]. The UPS code solves the PNS equations us- 
ing a fully conservative, finite-volume approach in a 
general nonorthogonal coordinate system. The UPS 
code has been extended to permit the computation of 
flowfields with strong upstream influences. In regions 
where strong upstream influences are present, the 
governing equations are solved using multiple sweeps. 
As a result of this approach, a complete flowfield can 
be computed more efficiently (in terms of computer 
time and storage) than with a standard N-S solver 
which marches the entire solution in time. Three it- 
erative PNS algorithms (IPNS, TIPNS, and FBIPNS) 
have been developed. The iterated PNS (IPNS) al- 
gorithm [14] can be applied to flows with moder- 
ate upstream influences and small streamwise sepa- 
rated regions. The time iterated PNS (TIPNS) algo- 
rithm [15] can be used to compute flows with strong 
upstream influences including large streamwise sep- 
arated regions. The forward-backward sweeping it- 
erative PNS (FBIPNS) algorithm [16] was recently 
developed to reduce the number of sweeps required 
for convergence. 

The majority of MHD codes that have been de- 
veloped combine the electromagnetodynamic equa- 
tions with the full Navier-Stokes equations result- 
ing in a complex system of eight scalar equations. 
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= V-(V-r)- V-U + E-J 


These codes can theoretically be used for any mag- 
netic Reynolds number which is defined as = 
VeVeV<x>L where a e is the electrical conductivity, \i t 
is the magnetic permeability, is the freestream 
velocity, and L is the reference length. However, 
it has been shown that as the magnetic Reynolds 
number is reduced, numerical difficulties are often 
encountered [4]. For many aerospace applications 
the electrical conductivity of the fluid is low and 
hence the magnetic Reynolds number is small. In 
these cases, it makes sense to use the low magnetic 
Reynolds number assumption and reduce the com- 
plexity of the governing equations. In this case, the 
MHD effects can be modeled with the introduction 
of source terms into the fluid flow equations. Several 
investigators [4, 8, 17-19] have developed N-S codes 
for the low magnetic Reynolds number regime where 
the induced magnetic field is negligible compared to 
the applied magnetic field. 

In the present study, a new PNS code (based on 
the UPS code) has been developed to compute MHD 
flows in the low magnetic Reynolds number regime. 
The MHD effects are modeled by introducing the ap- 
propriate source terms into the PNS equations. Up- 
stream elliptic effects can be accounted for by us- 
ing multiple streamwise sweeps with either the IPNS, 
TIPNS, or FBIPNS algorithms. Turbulence has been 
included by modifying the Baldwin-Lomax turbu- 
lence model [20] to account for MHD effects using 
the approach of Lykoudis [21]. The new code has 
been tested by computing both laminar and turbu- 
lent, supersonic MHD flows over a flat plate. Com- 
parisons have been made with the previous complete 
N-S computations of Dietiker and Hoffmann [18]. In 
addition, the new code has been used to compute the 
supersonic viscous flow inside a rectangular channel 
designed for MHD experiments [22]. 


( 3 ) 


Ohm’s law 

J = a e (E + V x B) (4) 

where V is the velocity vector, B is the magnetic 
field vector, E is the electric field vector, and J is the 
conduction current density. The flow is assumed to 
be either in chemical equilibrium or in a frozen state. 
The curve fits of Srinivasan et al. [23, 24] axe used 
for the thermodynamic and transport properties of 
equilibrium air. 

The governing equations are nondimensionalized 
using the following reference variables. 
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where the superscript * refers to the nondimensional 
quantities. In subsequent sections, the asterisks are 
dropped. 

If the flow variables are assumed to vary in only two 
dimensions ( x , y) while the velocity, magnetic, and 
electric fields have components in three dimensions 
(x, y, z), the governing equations can be written in 
the following flux vector form: 


Governing Equations 

The governing equations for a viscous MHD flow with 
a small magnetic Reynolds number are given by [18]: 
Continuity equation 

|e + v -(pv) = o a) 


au SFj 

dt + dx + dy 


dE v ,&F V 

-j^ + ^+Smbd 


( 6 ) 


where U is the vector of dependent variables and E* 
and are the inviscid flux vectors, and E v and F v 
are the viscous flux vectors. The source term Smhd 
contains all of the MHD effects. The flux vectors are 
given by 


Momentum equation 

^2+V- [pVV + p!] =V-f + JxB (2) 
Energy equation 

%^ + V-[(pe f +p)V] 


T 

u = [ p, pit, pv, pw, pe t ] ( 7 ) 
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where 


pe t = Ip ( u 2 + v 2 + w 2 ) + 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


and 7 can be determined from the curve fits of Srini- 
vasan et al. [23] for an equilibrium air flow or is equal 
to a constant (7) for a frozen or perfect gas flow. The 
nondimensional shear stresses and heat fluxes are de- 
fined in the usual manner [11]. 

The governing equations are transformed into com- 
putational space and written in a generalized coordi- 
nate system (£,77) as 

j Ut + E ? + F„ = ^p (13) 

where 

E = (Et — E„) + (jfj (F» - Ft,) 

F = ('J') (E« — E«) + (jj'} (F« — F«) (14) 


and J is the Jacobian of the transformation. 

The governing equations are parabolized by drop- 
ping the time derivative term and the streamwise 
direction (f) viscous flow terms in the flux vectors. 
Equation (13) can then be rewritten as 


where 

E 

F 


E< +F„ 


Smhd 

J 


(15) 


(7) E ‘ + (y) Fi 

(■y-) (E* ~~ E») + (-y) (F; - F(,) (16) 


The prime in the preceding equation indicates that 
the streamwise viscous flow terms have been dropped. 

For turbulent flows, the two-layer Baldwin-Lomax 
turbulence model [20] has been modified to account 
for MHD effects. Only the expression for turbulent 
viscosity in the inner layer is changed. This modifi- 
cation for MHD flows is due to Lykoudis [9,21]. 


Numerical Method 


The governing PNS equations with MHD source 
terms have been incorporated into NASA’s upwind 
PNS (UPS) code [13]. These equations can be solved 
very efficiently using a single sweep of the flowfield 
for many applications. For cases where upstream 
(elliptic) effects are important, the flowfield can be 
computed using multiple streamwise sweeps with ei- 
ther the IPNS [14], TIPNS [15], or FBIPNS [16] algo- 
rithms. This iterative process is continued until the 
solution is converged. 

For the iterative PNS (IPNS) method, the E vec- 
tor is split using the Vigneron parameter (cu) [25]. 
This parameter does not need to be changed for the 
present low magnetic Reynolds number formulation. 
In the previous high magnetic Reynolds number code 
[12] it was necessary to modify the Vigneron param- 
eter to account for MHD effects. After splitting, the 
E vector can be written as: 


where 


E* = — 
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E = E* 4- E p 


pu 

pu 2 + c op 
puv 
puw 

(pe t + p) u 

pv 
puv 

pv 2 + up 
pvw 

(pe t + p)v 


(17) 
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( 18 ) 


The streamwise derivative of E is then differenced 
using a forward difference for the “elliptic” portion 
(E*): 


(f ) . , = i? " e:) + (Ef+2 “ E? + i)] 

,+1 • ( 19 ) 
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where the subscript ( i 4- 1) denotes the spatial index 
(in the £ direction) where the solution is currently 
being computed. The vectors EJ +1 andEf +1 are then 
linearized in the following manner: 

E‘+i = E* + -U<) 

E?+i = E ?+(^) j (U i+ 1 -U 0 ( 20 ) 

The Jacobians can be represented by 


A* 

A p 


dE* 

au 

dE p 

dv 


( 21 ) 


After substituting the above linearizations into 
Eq. (19), the expression for the streamwise gradient 
of E becomes 



4-[(A*-A?)(U <+1 -U<) 

+ (Ef +2 - Ef) ] (22) 


The final discretized form of the fluid flow equa- 
tions with MHD source terms is obtained by substi- 
tuting Eq. (22) into Eq. (15) along with the linearized 
expression for the flux in the cross flow plane. The 
final expression becomes: 


Test Case 1: Supersonic laminar and 
turbulent flows over a flat plate with 
applied magnetic field 

In this test case, the supersonic, laminar and turbu- 
lent flow over a flat plate with an applied magnetic 
field is computed. This case corresponds to the flat 
plate case computed previously by Dietiker and Hoff- 
mann [18] using the full N-S equations. A strong 
magnetic field is applied normal to the flow as shown 
in Fig. 1. The dimensional flow parameters for this 
test case are: 


Moo 

= 

2.0 

P OO 

= 

1.076 x 10 4 5 N/m : 

Too 

= 

300 K 

R&oo 

= 

3.75 x 10® 

7 

= 

1.4 

L 


0.08 m 

<?e 

= 

800 mho/m 


The plate is assumed to be an adiabatic wall and a 
perfect gas flow is assumed. The magnetic Reynolds 
number (based on the length of the plate) is 0.056 and 
can be considered negligible when compared to one. 
The normal magnetic field component (B y ) ranges 
in value from 0.0 to 1.2 T. The magnitude of the 
magnetic field can be represented by the parameter 
m which is defined [18] by 





(Uj+i - U<) fc+1 


= RHS (23) 

where 

RHS = -i-[(Ef +2 )*-(E?)* +1 ] 

_ fd F\* +1 / SmhdN* 

\*i)i V j )i+i 

and the superscript k4*l denotes the current iteration 
(i.e. sweep) level. 


Numerical Results 

In order to investigate the utility and accuracy of 
the present PNS approach of solving MHD flowfields 
at low magnetic Reynolds numbers, a few basic test 
cases were computed. The supersonic viscous flow in 
these cases is altered by the presence of the magnetic 
and electric fields which are applied to the flow. 


PooU oo 


(24) 


and has units of (1/m). For B v = 1.2 T, m is equal 
to 1.33. 

A highly stretched grid consisting of 50 points in 
the normal direction was used to compute this case. 
The first point off the wall was located at 2 x 10 -7 
m. Initially, the flow was assumed laminar and sev- 
eral values of By ranging from 0.0 (no magnetic field) 
to 1.2 T were used. The velocity and temperature 
profiles at x = 0.06 m are shown in Figs. 2 and 3 for 
By = 0.0 T, 1.0 T, and 1.2 T. The velocity profiles 
are compared to the N-S results of Dietiker and Hoff- 
mann in Fig. 2 and show excellent agreement. The 
magnetic field generates a Lorentz force which acts in 
a direction opposite to the flow. Thus, the flow is de- 
celerated as the magnetic field is increased as seen in 
Fig. 2. For B y = 1.2 T the flow is slightly separated. 
The temperature profiles cannot be compared at this 
time since no temperature data is given in Ref. [18]. 

The turbulent flow over the flat plate was then 
computed using the modified Baldwin-Lomax turbu- 
lence model that accounts for MHD effects. The flow 
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was assumed laminar prior to the point (x = 0.04 m) 
where transition from laminar to turbulent flow was 
triggered. Again, several values of B y ranging from 
0.0 to 1.2 T were used in the computations. The tur- 
bulent velocity and temperature profiles at x = 0.06 
m axe shown in Figs. 4 and 5 for B y = 0.0, 1.0 T, 
and 1.2 T. The turbulent velocity profiles in Fig. 4 
are in good agreement with the results of Ref. [18]. 
The variation of skin friction coefficient is shown in 
Fig. 6. The present laminar /turbulent skin friction 
variations are compared with the results of Ref. [18] 
and show good agreement. The difference in results 
near the transition point may be due to the coarse 
grid and smoothing used in Ref. [18]. 

All of the present laminar computations were per- 
formed using a single sweep of the flowfield except for 
the separated flow case (B y = 1.2 T). For this case 
as well as for all the turbulent cases, multiple sweeps 
were used to account for upstream effects. 

Test Case 2: Supersonic viscous flow in 
a rectangular MHD accelerator 

In this test case, the supersonic flow in an experi- 
mental MHD channel is simulated. This facility is 
currently being built at NASA Ames Research Center 
by D. W. Bogdanoff, C. Park, and U. B. Mehta [22] 
to study critical technologies related to MHD bypass 
scramjet engines. The channel is about a half meter 
long and contains a nozzle section, a center section, 
and an accelerator section. The channel has a uni- 
form width of 2.03 cm. Magnetic and electric fields 
can be imposed upon the flow in the accelerator sec- 
tion. A schematic of the MHD accelerator section is 
shown in Fig. 7. 

This test case was previously computed by 
R. W. MacCormack [10] using the full N-S equations 
coupled with the electromagneto dynamic equations. 
The electrical conductivity in his calculations was set 
at 1.0 x 10 5 mho/m resulting in a very large magnetic 
Reynolds number. In the present study, the calcu- 
lations are performed in the low magnetic Reynolds 
number regime using a realistic value of electrical con- 
ductivity. The flow is computed in two dimensions, 
but later will be extended to three dimensions. Be- 
cause of flow symmetry, only half of the channel is 
computed in the 2-D calculations. 

The flow in the nozzle section and the center sec- 
tion was computed using a combination of the OVER- 
FLOW code [26] and the present PNS code (without 
MHD effects). The initial conditions for the nozzle 
(flow at rest) were: 

PO = 8.0 x 10 s N/m 2 


To = 7500 K 

The laminar flow was assumed to be in chemical equi- 
librium. The computed flowfield at the end of the 
center section was then used as the starting solu- 
tion for the flow calculation of the accelerator section. 
The MHD parameters used in the accelerator section 
were: 


<y* 

= 50 mho/m 

By 

= 1.5 T 

E z 

= -Ku c B y 

Be m 

= 0.05 


where the load factor ( K ) ranged in values from 0.0 to 
1.4, and the centerline velocity (it c ) at the beginning 
of the accelerator section had a value of 3162 m/s. 

The velocity profiles at the end of the accelerator 
section are shown in Fig. 8 for different load factors. 
The velocity profile with no electric or magnetic fields 
is denoted by K = 0. The increase in the centerline 
velocity with distance (a;) for various load factors is 
shown in Fig. 9. The centerline velocity increases 
by about 30% with a load factor of 1.4. It should 
be noted that the flow decelerates because of friction 
when no electric or magnetic fields are applied. 

Concluding Remarks 

In this study, a new parabolized Navier-Stokes al- 
gorithm has been developed to efficiently compute 
MHD flows in the low magnetic Reynolds number 
regime. The new algorithm has been used to com- 
pute both laminar and turbulent, supersonic, MHD 
flows over flat plates and in a rectangular accelera- 
tor section. Although only limited results have been 
obtained thus far, it can be seen that the present 
approach is quite promising. Computations of other 
test cases axe currently underway in order to validate 
the current method. 
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Figure 7: Schematic of MHD accelerator section 
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Figure 9: Centerline velocity for various load factors 
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